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Rotation of non-circular hydraulic jumps is a recent experimental observation that lacks a theory based on 
first principles. Here we furnish a basic theory of this phenomenon founded on the shallow-water model of the 
circular hydraulic jump. The breaking of the axial symmetry morphs the circular jump into a polygonal state. 
Variations on this state rotate the polygon in the azimuthal direction. The dependence of the rotational frequency 
on the flow rate and on the number of polygon vertices agrees with known experimental results. We also predict 
how the rotational frequency varies with viscosity. Finally, we establish a correspondence between the rotating 
polygonal structure and the Hawking effect in an analogue white hole. The rotational frequency of the polygons 
affords a direct estimate of the frequency of the thermal Hawking radiation. 
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In low-dimensional flows, the hydraulic jump is associated 
with an abrupt increase in the height of a flowing liquid jl]]. 
Along a standing circular locus, this feature is known as the 
two-dimensional circular hydraulic jump Bll. Its geome¬ 
try maintains a symmetry about a reference axis that passes 
normally through the point of the origin of the flow on a flat 
plane. To develop a cogent mathematical theory, the flow is 
viewed as a shallow layer of liquid E confined to a horizontal 
plane, diverging radially outwards from a source point. There¬ 
after, applying the boundary-layer approximation in the shal¬ 
low flow, viscosity is invoked |^,|4{], along with all other com¬ 
plications of the nonlinear Navier-Stokes equation. A com¬ 
pelling evidence in favour of viscosity comes from an experi¬ 
mentally verified formula of the jump radius, scaled in terms 
of viscosity ES. Various physical interpretations of this 
scaling law have been forwarded ESl- 

Circular hydraulic jumps appear in Type-I and Type-II 
states |@], of which the latter is formed by increasing the 
height of the flow in the post-jump region. Type-II jumps have 
a wider jump region with a surface eddy all along its circular 
rim like a floating torus @]. When water is replaced with a liq¬ 
uid of much greater viscosity, the Type-II state spontaneously 
breaks the axial symmetry of the circular state, resulting in 
a front that has a polygonal geometry ffl. This transition 
has visibly temporal features because a linear instability arises 
from the initially stable circular state, and leads to the for¬ 
mation of a non-axisymmetric polygonal structure, whose de¬ 
pendence on the azimuthal angle is periodic E- Within the 
widened jump region, the transport of liquid also becomes az¬ 
imuthal, although prior to it, in the pre-jump region, the flow 
is radial E- Clearly, the breaking of axial symmetry is lo¬ 
calized only at the jump radius or thereabouts. All of these 
features of the polygonal jump have been known well for two 
decades, but in addition, as reported very recently, the polyg¬ 
onal jump undergoes a rotation 0. This new phenomenon 
has been named the “rotational hydraulic jump” [10] 12411 . 


Thus far we have provided a brief narrative of the role 
that viscosity plays in the multifarious properties of the two- 
dimensional hydraulic jump. Now we consider the jump from 
a different perspective [3,11]- 14], according to which the po¬ 
sition of the standing jump is a boundary where the steady ra¬ 
dial velocity, vo (r), equals the local speed of long-wavelength 
surface gravity waves, y/gho, with ho(r) being the steady lo¬ 
cal height of the flow layer. This boundary is like a standing 
horizon, segregating the supercritical and the subcritical re¬ 
gions of the flow, with criticality referring to the condition 
under which the speed of the bulk flow matches the speed 
of surface gravity waves E- This entire point of view lays 
strong emphasis on the advective and pressure terms in the 
Navier-Stokes equation, to the neglect of the viscous term. 
The radius of the horizon (which is also the jump), rj, is de¬ 
fined by the critical condition, Vg(rj) = gho(rj), and sets a 
spatial limit for the transmission of information. As the equa¬ 
torial flow proceeds from its point of origin, its radial velocity, 
vo > y/gho, but viscosity and the radial geometry slow down 
the flow. When the critical condition is met, both the jump 
and the horizon occur simultaneously E- In the supercritical 
part of the jump, where vo > y/gho , gravity waves (as car¬ 
riers of information) cannot travel upstream against the bulk 
flow, and hence every point in the supercritical region remains 
uninformed about the fate of the flow downstream. This state 
of affairs prevails everywhere within the jump, and so it acts 
as an impenetrable barrier against the percolation of any in¬ 
formation from the outside to the inside. The horizon implied 
by the critical condition has been amply demonstrated in, 
but the horizon by itself is inadequate to explain why a jump 
should coincide with it, a point that has been qualitatively ad¬ 
dressed and appreciated for long E- Concisely stated, the 
horizon is a necessary condition for the jump but not a suffi¬ 
cient one, and so without reference to anything regarding the 
jump, the horizon is just an analogue of a white hole. In this 
study we show that the formation of a polygonal jump and 
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its observed rotation are the natural outcomes of break¬ 
ing the axial symmetry of the circular horizon of the analogue 
white hole. The breaking of the axial symmetry, localized 
near the horizon, is due to the analogue surface gravity. The 
rotation has a persuasive similarity with the Hawking effect, 
which is also due to the analogue surface gravity. 

The mathematical description of the two-dimensional flow 
is most succinctly framed in the cylindrical coordinate sys¬ 
tem, (r, <f>, z), and by tailoring the Navier-Stokes equation 
accordingly [1]]. Our analysis starts with the steady height- 
integrated Navier-Stokes equation of a shallow-water radial 
base flow J2l, 


vo- 
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dr 


dh 0 v 0 

g = -v 


dr 




(1) 


and the steady height-integrated equation of continuity. 
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in its differential form. The integral form of Eq. (O is rvo ho = 
Q/2ir, where Q, a constant, is the steady volumetric flow rate. 
The subscript “0” in the foregoing equations stands for steady 
radially varying quantities, of which the velocity of the flow, 
vo (r), has been obtained in the shallow-water theory by ver¬ 
tically averaging the radial component of the velocity across 
the height of the flow. The boundary conditions used for the 
averaging are that velocities vanish at z = 0 (the no-slip con¬ 
dition), and vertical gradients of velocities vanish on the free 
surface of the flow (the no-stress condition) 62, 14, 3. These 
boundary conditions are applied under the standard assump¬ 
tion that while the vertical velocity is much small compared 
with the radial velocity, the vertical variation of the radial ve¬ 
locity (through the shallow layer of water) is much greater 
than its radial variation j3|. In this manner all dependence on 
the z-coordinate is averaged out. 

Now, about the steady radial background flow, as implied 
by Eqs. ([TJ and (0, we develop a time-dependent azimuthal 
perturbation scheme prescribed by h(t,r,cp) = hp(r) + 
h'(t, r, (j>), v r (t, r, <p ) = v 0 (r) + v’ r (t, r, <j>) and v${t, r, <p) = 
0 + r. <fi). All primed quantities are time-dependent per¬ 
turbations in both the radial and azimuthal coordinates, with 
v' r and v' 0 being perturbations on the radial and azimuthal ve¬ 
locity components, respectively. We have designed the az¬ 
imuthal flow to be entirely a perturbative effect without any 
presence in the steady background flow. Our next step is 
to define a new variable, / = rhv r , whose steady value, 
fo = rhovo, is a constant, as we can see from Eq. 0. Within 
a multiplicative factor, owed to the two-dimensional geometry 
of the system, fo gives the conserved volumetric flow rate un¬ 
der steady conditions. So f is a perturbation on this constant 
background volumetric flow rate. Linearizing in /' under the 
formula, f = fo + /', gives us /' = r (hov' r +Voh'). A 
similar linearization of the general time-dependent continu¬ 
ity equation, bearing both radial and azimuthal variations jTj|, 


leads to 

dh! 1 df ho dv^ 

dt r dr r dp 


(3) 


Likewise, from the radial and azimuthal components of the 
Navier-Stokes equation, as expressed in cylindrical coordi¬ 
nates i, we derive 


dv f d 


ah' 
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and 
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r dp 
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respectively. We stress here that in extracting Eqs. 0 and 0, 
we have ignored all product terms of viscosity and the primed 
quantities. Therefore, in our treatment, viscosity makes its 
mark through the zero-order stationary quantities, vo and h q, 
as Eqs. 0 and 0 indicate. In all first-order equations in¬ 
volving the primed quantities, Eqs. 0, 0 and 0, viscosity 
does not appear explicitly, and exerts its influence implicitly 
through the zero-order coefficients. This is a satisfactory ap¬ 
proximation, to the extent that our principal concern is the 
effect of azimuthal variations on the axially symmetric radial 
flow, and the breaking of the axial symmetry therefrom. 

We use Eq. 0 to substitute dh'/dt in the first-order partial 
time derivative of /', whereby we also obtain 


dv' r __ vp df vl df h 0 vl dv 0 
dt f 0 dt f 0 dr f 0 dp ' 

Collectively, Eqs. 0 and 0 present a closed set of condi¬ 
tions that express h' and v' r as a linear combination of f and 
Vc/,. At this stage we require two independent mathematical 
conditions on which we can impose Eqs. 0 and 0. Such 
conditions are readily supplied by the perturbations of the ra¬ 
dial and azimuthal dynamics, as shown in Eqs. 0 and 0, 
respectively. We take the second-order partial time derivative 
of these two coupled equations, and on them we apply the con¬ 
ditions provided by Eqs. 0 and 0, along with the second- 
order partial time derivative of Eq. 0. This long algebraic ex¬ 
ercise ultimately delivers two equations that are second-order 
in time. They are 
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derived from the perturbation of the radial dynamics, as given 
by Eq. 0, and 


d 2 v <j) vo d f dv^\ ghod 2 v lj) _ g d fdf\ 

dt 2 r dr \ dt J r 2 dip 2 r 2 dip \ dr J ’ 












3 


derived likewise from the perturbation of the azimuthal dy¬ 
namics, as given by Eq. (0). The two foregoing equations, 
linearized and coupled, form a set of wave equations in f and 
A familiar wave equation in /' only, pertaining just to the 
radial dynamics, is obtained when all the (/(-derivatives on the 
right hand side of Eq. CD are made to vanish jsj]. Thereafter, 
the expression on the left hand side of Eq. CD is rendered com¬ 
pactly as d a (f a/3 dpf) = 0, in which the Greek indices run 
from 0 to 1, with 0 standing for t and 1 standing for r. This 
simplification establishes an acoustic metric and an acoustic 
horizon in the physical problem of the two-dimensional hy¬ 
draulic jump, the details of which have been discussed in a 
previous work |3|]. Of course, we must remember that this 
reasoning is valid only for an inviscid fluid in a potential flow. 
In our present approach, the steady background flow is af¬ 
fected by viscosity, while the first-order perturbation has an 
azimuthal component, which we consider to be a major ad¬ 
vancement. It cannot be ignored at the jump front, as far as 
the formation of polygonal structures is concerned. 

We anticipate solutions that are separable in f, r and 
(/>, for the two linearly coupled wave equations, given by 
Eqs. © and 0. In keeping with this stipulation, we set down 
</>) = Aexp[— iuut + i<d(r) + imp] and ^(f, r, (/>) = 
Bip(r) exp [—iu)t + i<d(r) + imp] = Bf'ip{r)/A , in which 
A and B are constants. The factor of captures the az¬ 
imuthal variation. Along the radial direction, /' has a slow 
variation, except near the acoustic horizon. The horizon is sig¬ 
nificant because it is also the position of the hydraulic jump, 
where the axial symmetry of the circular jump front is bro¬ 
ken, a polygonal structure emerges and the flow acquires an 
azimuthal component. Accordingly, v<f, is most conspicuous 
about the jump radius and decays sharply away from it, a fea¬ 
ture that is captured by a strong local prominence for 'p(r) at 
r = rj. This radial position also allows us to exploit the con¬ 
dition of the acoustic horizon, Vq = gho, greatly simplifying 
our analysis. And so Eqs. 0 and (0 yield two characteristic 
equations as 
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respectively. We emphasize that vo(rj), ho(rj) and their 
derivatives in Eqs. 0 and (flOt carry their values only at 
rj. The consistency of Bp/A on the left hand sides of 
Eqs. 0 and (flOb . which are both quadratic in uj, lets us elimi¬ 
nate A/B, resulting in a single quartic equation in cu, going as 
(uj 3 +T 2to 2 +Tiu; + To)u; = 0. One root of the quartic equa¬ 
tion is uj = 0 , which leaves behind a residual cubic equation. 


Marginal stability (signalling the onset of a possible instabil¬ 
ity) extracts another root of w = 0 from the cubic equation, 
something that is possible only when T 0 = 0 , and which in 
turn leads to a critical value of m = m c . From the vanishing 
of the real part of Tq we get 
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and similarly the vanishing of the imaginary part of Tq gives 
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which determines d0/dr at the horizon. 

Having obtained m c from the condition of marginal stabil¬ 
ity, we now cause a slight perturbation, m = m c + Am (with 
A m/m c < 1 ), to effect a very small change in uj from u> = 0 . 
Smallness of the value of uj allows us to ignore its higher or¬ 
ders in the cubic equation and retain only Tiw + To ~ 0 . 
Hence, — u> ~ To/T!. On linearizing in Am, we get 


r 0 = 


2m c gh 0 v 0 


d0 . 2 dt'o \ 
dr vq dr J 


Am, 


(13) 


while we write Id = P + iQ, in which, aided by some sim¬ 
plifications due to Eqs. £jT]> and dm 


P = v 2 


d0 

dr 


d (gho) /dr 


Q = - 


d0 

dr 


d{vl)/dr 

d ( 2v 2 -gho) + ^ 
r j 


dr 


1 + 


rj dip\ 
ip dr J 


For convenience we write — u> = 7 Am, from which we get 


2m c v 0 gh 0 
1 ~ r 2 ( P 2 + Q 2 ) X 

pd0 _ 2Q_ dro\ _ . / d0 | 2P dv 0 
dr vo dr J \ dr vq dr 

whose form is best read as 7 = 5ft(7) + *3( 7 ). Now, the 
phase of the wave solution is exp[— iu>t + i@(r) + im<j>], from 
which we extract only the time-dependent part and recast it 
as exp[i 3 ?( 7 )Amf] x exp[— 3 ( 7 )Am£]. The conclusion we 
draw is that 3 ^( 7 ) causes the wave to travel along the <p coor¬ 
dinate (the azimuthal direction), and 3 ( 7 ), depending on its 
sign, causes either a growth or a decay in the amplitude of 
the travelling wave. To examine the latter feature, we look at 
Eqs. (ITTb and (fl2] >. both of which give m 2 as a perfect square. 
As such, m c will have two roots of the same value but op¬ 
posite signs. Since 3 ( 7 ) depends on m c , as Eq. (ITTb shows, 
it will similarly carry both signs, with the signs of all other 
quantities in Eq. (ITTb arguably being fixed. If 3( 7 ) > 0, then 
stability will be achieved only by Am > 0, i.e. if m > m c . 
This is to say that a polygon, so formed, will be stable only if 
the number of its vertices is above a threshold given by m c . 
The opposite argument applies if 3 ( 7 ) < 0, because in this 
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case stability is ensured by m < m c , with there being an up¬ 
per limit to the vertices of a stable polygon. 

Recent experiments by Teymourtash and Mokhlesi E3i 
have shown unstable polygons to undergo a rotational be¬ 
haviour. We reproduce this feature theoretically with the help 
of Jf( 7 ), in which the sign of m c controls the clockwise or the 
anticlockwise direction of the rotation. The angular frequency 
of the rotation is fl ro t = |5R(7 )At7i|. In an unstable situation 
where m < m c , small values of m yield high values of | Am|, 
and so fl rot oc | m c — m\. This linear decay of the angular 
frequency with increasing number of vertices matches exper¬ 
imental results lIToh . To find how O rot depends on the flow 
rate, Q, we first note from Eq. (IT4l) that Q rot ~ vo{rj)/rj, 
in which rj ~ Q b / 8 i/~ 3 ^ 8 g~ 1 / 8 , a well-known scaling re¬ 
sult 121 that is derived by equating the dynamic time scale 
of the steady radial flow with the time scale of viscous dis¬ 
sipation 0. The steady background quantities depend on 
viscosity, and at the jump, where tto(rj) = \/gho(rj), the 
flow height, ho(rj ), is scaled by combining the aforemen¬ 
tioned scaling of rj with the integral solution of Eq. ©. 
This results in /io(rj) ~ Q 1 / 4 !/ 1 / 4 ^- 1 / 4 with which we 
get the scale, Q rot ~ Q~ x l 2 v x l 2 g 3 l 2 . Evidently, Q rot de¬ 
creases with increasing flow rate, Q, something that has been 
observed by Teymourtash and Mokhlesi m in their experi¬ 
ments. Our theoretical treatment is, therefore, well in accord 
with two experimental findings of Teymourtash and Mokhlesi 
Ei. namely, the two ways for fl rot to decay — with increas¬ 
ing number of polygon vertices and with increasing flow rate. 
Beyond these two established facts, we make a prediction, 
based on O rot oc v 1 ! 2 , that the angular velocity of the rotating 
polygons will increase with increasing kinematic viscosity. 

Our most crucial claim is that the breaking of the axial 
symmetry has a connection with the Hawking radiation in an 
acoustic white hole. If the right hand side of Eq. 0 were 
to vanish, the left hand side will bring forth the symmet¬ 
ric metric of an analogue white hole yl 111 11411 . a point of 
view that is applicable to a steady radial outflow. The in¬ 
dependent terms on the right hand side of Eq. 0 break the 
axial symmetry of the steady radial flow, and at the hori¬ 
zon, where Vq = ghg, the flow becomes azimuthal. The 
analogue “surface gravity” at the horizon EM] is given as 
g s = (—1/2) x [d(t?g — gho)/dr], which, taken together with 
Eq. ( ITTb . means m 2 ~ [gs/fo^J)]^rot- Further, in terms of 
g s , the Hawking temperature, Th, can be set down as j 1911 
^bTh = Ml h, in which Qh = g s /[27rr;o(rj)] is a charac¬ 
teristic frequency of the thermal Hawking radiation. Clearly, 
H rot ~ Qn, whose startling implication is that the frequency 
with which the polygon rotates is a practicably measurable 
manifestation of the Hawking effect in an analogue white 


hole, a statement that is in conformity with a very recent ob¬ 
servation of Hawking radiation emanating from an analogue 
black hole in an atomic Bose-Einstein condensate 0. The 
corresponding Hawking temperature specifies a thermal scale 
for any energy involved in the Hawking process [21], and at 
high frequencies, w fin, the tunnelling amplitude of the 
Hawking radiation I 22 I 1 assumes the recognizable form of the 
Boltzmann factor, exp[— (Mc)/(£;bTh)]. 
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